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We investigate the memory effect in a simple model for glassy relaxation, a trap model with a 
Gaussian density of states. In this model thermal equilibrium is reached at all finite temperatures 
and therefore we can consider jumps from low to high temperatures in addition to the quenches 
usually considered in aging studies. We show that the evolution of the energy following the Kovacs- 
protocol can approximately he expressed as a difference of two monotonously decaying functions 
and thus show the existence of a so-called Kovacs hump whenever these functions are not single 
exponentials. It is well established that the Kovacs effect also occurs in the linear response regime 
and we show that most of the gross features do not change dramatically when large temperature 
jumps are considered. However, there is one distinguishing feature that only exists beyond the 
linear regime which we discuss in detail. For the memory experiment with 'inverted' temperatures, 
i.e. jumping up and then down again, we find a very similar behavior apart from an opposite sign 
of the hump. 



I. Introduction 

The out-of-equilibrium dynamics of glassy systems has been studied for a long time, for a review 
see[l]. Many investigations have been devoted to the study of the violations of the fluctuation- 
dissipation theorem, cf. ref. [2]. In particular, the question regarding a meaningful definition of 
an effective temperature still remains to be answered as there are conflicting results, for a recent 
investigation see|3j and the literature therein. Also the aging behavior of time-correlation functions 
and in particular the rotational and translational dynamics has been investigated with the result 
that the aging properties of glass-forming liquids are determined by the primary relaxation, cf. 
ref. HE]. 

Another effect that has been analyzed in much detail both theoretically and experimentally in 
the past decade is the so-called Kovacs effect or memory effect [HI [7j. In the original version of 
the experimental protocol one equilibrates a glass-forming liquid (low molecular or polymeric) at 
a given temperature Tq. Afterwards, a temperature jump to a lower temperature Ti is performed, 
followed by an increase in temperature to a final temperature T2 after a time ti and the volume 
of the sample is monitored as a function of time. The crucial point in the experiment is that 
the time ti for the up-jump Ti — )■ T2 is chosen such that Vri(ti) coincides with the equilibrium 
volume at T2, VTi(ti) = V^'^. Instead of staying constant at Vj^', the volume first increases and 
then decreases again as a function of time and V^^ is only reached in the long-time limit. The 
observation of this so-called Kovacs hump usually is taken as an indication for the fact that the 
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system cannot be described fully in terms of time-independent thermodynamical quantities. The 
relaxation appears to depend sensitively on the history of the sample. This memory effect has not 
only been observed in the evolution of thermodynamic variables but also shows up in the dielectric 
constant in both, polymer glasses[8] and a molecular liquid|t9j. Furthermore, its occurence is not 
restricted to the slow time scale of the primary relaxation but it has also been found in the 
region of the low-frequency vibrations of a glass[10j. The Kovacs effect and its interpretation in 
terms of simple models has been described in a review-like paper by Bertin et al. [TT] . Among the 
models investigated is the so-called trap model with an exponential density of states (DOS) [12]. 
Also computer simulations of model glass-forming liquids [13] and the Ising spin glas[Tl] show 
the Kovacs effect as well as a number of models showing glassy relaxation, such as the p-spin 
modellblSj, facilitated or constrained spin models [151 [171 HH] and also the two-level model[in]. The 
effect has furthermore been observed in other complex systems like granular materials [20] and also 
in a model for critical dynamics |21j. 

If one chooses to perform the experiment with electrical fields as a variable instead of the 
temperature, one can perform so-called crossover experiments in the linear response regime and 
it was shown that the hump has its origin in the fact that one subtracts two relaxation functions 
which are evaluated at different times [22]. This means, one has to build +t) — 62^{t) with 

denoting the dielectric (step) response function and the 6k are the relative field intensities. 
Thus, in the linear regime, the Kovacs effect should be observable in all systems exhibiting non- 
exponential relaxation. Of course, this finding has its origin solely in the superposition principle 
of linear response theory and here the Kovacs hump is not at all related to any thermal history. 
A discussion of the behavior of the Kovacs hump in linear response theory for Markov processes 
has been presented recently [23]. 

In the present paper we study the Kovacs effect in the trap model with a Gaussian DOS. In 
contrast to the model with an exponential DOS, the system can reach equilibrium in the long 
time limit for all temperatures [24] . Additionally, this model has been shown to qualitatively 
reproduce some relevant features of the energy landscape of simulated glassy systems [251 126] . In 
particular, the relaxation time for a given energy has been shown to coincide with the escape rate 
from a meta-basin up to a factor of two. Furthermore, after a temperature jump the temporal 
evolution of the distribution of trap energies behaves similar to what has been observed in computer 
simulation studies of the distribution of inherent structure or meta-basin energies during aging 
of a model glass-forming liquid[27l I2B]- After a quench, the width of the distribution decreases 
and then increases again as a function of the time that has elapsed and for a jump from low 
to high temperatures one finds a bimodal distribution of meta-basin energies as predicted by the 
Gaussian trap model[21]. Below, we will briefly recall the non-exponential relaxation of the energy 
after both, a quench from a high temperature and an up-jump from a low temperature. Also for 
the Kovacs protocol we will in addition to the standard procedure treat the situation in which 
the system is in equilibrium at low temperatures initially. Then a jump to a higher temperature 
and afterwards back to some intermediate temperature is performed and the energy is monitored 
as in the original protocol. Furthermore, we will consider the Kovacs experiment for very small 
temperature differences in the linear response regime. 
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II. Non-exponential energy relaxation in the Gaussian trap 
model 



Before we discuss the Kovacs effect in the Gaussian trap model in detail, we briefly review the 
aging-properties of the model |29j and the relaxation after a temperature-jump [25] . 

The trap model is defined by the following master equation for the conditional probability 
(Green's function) of finding the system in a trap with trap energy e at time t + to, given it 
occupied the trap with energy eo at time t = Iq {GT{€,t + tokoi^o) = GT{€,t\eo,0) = GT{€,t\eo) 
due to time-translational invariance) [121 1211 [281 129] : 

Grie, t\eo) = -/€T(e)G'r(e, t\eo) + p{e) J dt\T{t)GT{t\ t\t^) (1) 

The subscript T is used to emphasize that the Green's function is defined at temperature T = 1/(5. 
In eq.([T]), the escape rate is given by 

Kr(e) = ^^ooe^' (2) 

with the attempt rate /too- We will measure time in units of /too in all following calculations. 
Furthermore, we solely consider the model with a Gaussian DOS 

p{e) = ^e-^'/^^'^'^ (3) 
y Z-na 



with 0" = 1 and throughout the paper we measure temperature in units of a. From eq.Q, 
the equihbrium populations at a given temperature T are found to be Gaussian P7?(e) = 
limt_s.oo GT(e, t|eo) = 75^^"''''"''^^ ^''^'^ ■* with It = —jSa'^. The solution of eq.(l| allows the com- 
putions of all quantities of interest. For instance, the population of the trap with trap energy e is 
given by 

Prie, ^) = J deoGrie, t - to|eo)p(eo, ^o) (4) 

where the integral Jde is to be read as a short hand notation for J^^de and p(eo,to) denotes the 
population at the initial time to. In all calculations performed in the present paper we assume that 
an instantaneous temperature jump from Tq to T took place immediately before to- This means 
that the initial population at the destination temperature is given by the equilibrium population 
at the reference temperature Tq, 

p(eo,io) =PTo(eo)- (5) 

Eq. Q allows the computation of general expectation values, for instance moments of the energy- 
distribution: 

E^it)= /rfee>T(e,t)- (6) 



The relaxation behavior of the energy has been reviewed briefly in ref.[28]. The mean value Exit) 
changes monotonously from the equilibrium value at the initial temperature Tq, E{t = 0) = E^ = 
— /3oO"^ towards the equilibrium value at the destination temperature T, E{oo) = E^ = —f3a'^. 
Also the change of the width of the distribution of the populations of the traps for intermediate 
times and in particular the intermittant two-peak structure after a temperature up-jump has been 
discussed. For our present purpose these latter observations are more important because the 
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change in the width of the populations pri^, t) gives rise to varying degrees of non-exponentiahty 
in the energy relaxation. This feature becomes apparent when the normalized relaxation function 



= (Erit) - Et{oo))/{Et{0) - ^t(oo)) 
with Et{0) = — /SqC"^ and Et{oo) = —fia"^ is recorded, cf. Figjlj One can see that the relaxation 




t 

Figure 1: Relaxation function versus time (in units of Kqo) for different temperatures (in 

units of 0"). Upper panel: jump from Tq = 0.5 to T = 0.45,0.4,0.35,0.3,0.25 from left to right. 
Lower panel: jump from Tq = 0.2 to T = 0.45, 0.4, 0.35, 0.3, 0.25 from left to right. 

time increases with decreasing temperature T independent of the initial temperature Tq and the 
non-exponentiality of the decay becomes more pronounced. For the jump from a low temperature, 
this behavior might be counter-intuitive on first glance. However, it just reflects the relaxation 
behavior at the destination temperature and this relaxation is faster the higher the temperature 
is. 

The relaxation times te{T,Tq) defined as the 1/e-decay times of $_B(t) for a variety of initial 
temperatures Tq as a function of the final temperature T are shown in Fig|2} It is apparent that for 
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Figure 2: te(T,To) for Tq = oo, Tq = 1, Tq = 0.3, Tq = 0.25 and Tq = 0.2 from bottom to top as 
a function of the inverse destination temperature /3 = T~^. Additionally shown is the equilibrium 
value Teq (dashed line) and the stretching parameter Xeq (inset) obtained from a fit of Ceq(t) to a 
KWW-function (exp [-(t/xeq)^^'']). 

a given temperature T the value of the relaxation time strongly depends on the initial temperature 
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and that te increases with decreasing To as is also observed from FigjTj This is just the same 
behavior as is observed if te is considered as a function of T. In addition to the relaxation times 
of we show the relaxation time Tcq of the normalized energy autocorrelation function Ceq(t) 

geq(t) = ~ With (e(t)e(0)) = |rfe|rfeoeeoGT(e,t|eo)p?>(eo). (7) 

The relaxation time is obtained from a fit to a KWW-function, Ceq(t) ~ exp [— (t/rcq)^'"^]. The 
energy autocorrelation function Ceq(t) is of relevance because it determines the energy relaxation 
in the linear response regime, which is explicitly given by ET{t) = + {^t[, ~ -^t*) C'eq(i)- As it 
is always the lower temperature that determines the value of te{T,Tq) it is not astonishing that 
Teq has values similar to the corresponding te- The stretching shows the same qualitative behavior 
as also observed for other observables in the framework of the Gaussian trap model |12l [29] . 

Instead of further discussing the relaxation properties of the Gaussian trap model after fast 
temperature jumps, we refer the interested reader to refs.p^ l27t [28] where a detailed discussion 
of the temporal evolution of the trap-populations pT{e,t) can be found. We close the discussion 
of the energy relaxation in summarizing that the energy relaxes towards equilibrium in a highly 
non-exponential way with both, the relaxation time and the stretching strongly depending on 
temperature. These features will be essential for our understanding of the origin of the so-called 
Kovacs hump. 

III. Kovacs-EfFect 

We will now discuss the memory effect for two situations, first starting from a high temperature 
and afterwards we consider the 'inverse' situation starting from a low temperature. In both cases 
we will consider the usual situation in which the temperature differences are not restricted to be 
small. The relevant questions that we will consider are whether the results for the Gaussian trap 
model are compatible with the features observed in more complex models and if the two-peak 
structure observed in the distributions of trap energies after an up-jump has some impact on the 
results of the 'inverse' Kovacs experiment. Finally, we will discuss the effect of the experiment in 
the linear response regime in order to see if there are intrinsic nonlinear effects that give rise to 
distinctive features for large temperature jumps. In all cases we choose the energy as the relevant 
observable. 

1. 'Standard' Kovacs experiment 

We start from a high initial temperature and follow the usual protocol, cf. Figj3| 

• quench: Tq — )■ Ti (Ti < T2 < Tq) 

• after a time ti chosen such that = : jump to T2. 

• measurement of the energy ET^it) (time t = at the beginning of the measurement). 
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Figure 3: Sketch of the experimental protocol for the Kovacs effect: In the beginning, a quench 
from To to Ti is performed. Then the system is kept at Ti until the energy has reached the 
equilibrium value at a temperature T2 (Ti < T2 < Tq). At this time ti, the temperature is 
instantaneously changed to T2 and the temporal evolution of the energy is monitored. 



From eqns. and (§) one finds for the temporal evolution of the energy 

ET^ih + t) = Jde-e-pT2{e,ti;t) = Jdejde ■e-GT2{(^,t\e)pT^{e ,ti] 



(8) 



where the populations at Ti follow from pxii^', h) = Jde^GT^^e', ti|eo)p2?(eo)- (Note that according 
to our notation defined in eq.Q one has PTii^-, h', t) = PXii^, h + ^)-) 

The definition of ti, ET^{ti) = E^, makes clear that Ex^iti + 1) can conveniently be written 

as 

ETAti + t) = E'T,+^E{t) (9) 

thus defining the quantity A£'(t), the so-called Kovacs hump. For this, one has by definition 
A£^(0) = and additionally lS.E{t — t- 00) = 0. For a more detailed discussion of the above 
definition see ref.|ll]. 

We start the discussion with the determination of the time ti defined by ET^{ti) = E^. The 
dependence on the temperature-difference T2 — Ti of ti is shown in Fig|4] starting from Tq = 00 
and To = 1 for different T2. The most important feature that becomes evident immediately is the 
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Figure 4: ti for T2 = 0.5, 0.6, 0.7 following a quench from To = 00 to various Ti G (0.2 
lines) and the same for a quench from To = 1 (dashed lines). 
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non-monotonic behavior of ti as a function of T2 — Ti which to the best of our knowledge has not 
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been observed previously in theoretical or experimental investigation of the Kovacs effect. This 
non-monotonic behavior might be rationalized if one compares ti to the times required to reach 
thermal equilibrium at Ti, Tcq. For the temperature range considered, Teq changes by more than 
four orders of magnitude, cf. Fig|2j Therefore, we interprete the non-monotonic behavior in the 
following way: 

Ti <^ T2: ti is long, because Teq(Ti) is long. Thus, it takes a long time to reach Ex^iti) = E^. 

Ti — > T2: Here, and E^ take on very similar values. Accordingly, Ex^it) has to relax 
almost to equilibrium in order to reach E^. In the limit Ti — )■ T2 ti diverges. 

Of course, the temperature at which the minimum value of ti is observed depends on the stretching 
of the relaxation function and also on the absolute relaxation strength {E^ — E^ in our case). 
Therefore, it is not easy to predict how large a temperature range has to be considered in a real 
experiment. We speculate that the reason why this anomalous behavior has not been observed 
up to now is that one has to use quite large temperature differences. This, however, should not 
provide a serious problem in possible experimental realizations or computer simulations. 

Next, let us consider the Kovacs-hump AE{t) as defined in eq.([9]). This quantity evolves as 
expected and as has been found in many earlier investigations. For Tq = 1 — )■ T2 = 0.6 it is plotted 
in FigjSj It is evident that for increasing Ti the hump shifts towards longer times and becomes 
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Figure 5: 'Kovacs-hump' vs. t; Tq = 1 — )■ T2 = 0.6 and various Ti, T\ = 
0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.575 from top to bottom. 



less intense. Therefore, the maximum 

AE^= AE(t)|_ (10) 

increases as a function of the temperature difference (T2 — Ti) and the time of the occurence of 
the maximum, 

Tk : AE{t = Tk) = AEk (11) 

decreases. 

Figj6] shows AEk and tk plotted versus (T2 — Ti). It is obvious how AE^ increases and 
decreases almost power-law like as a function of the temperature difference as also observed for 
other mo dels [T^. 
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Figure 6: AEk and as a function of (T2 — Ti); the lines are power laws AE^ — a ■ (T2 — Ti) 
and Tk ^ V(^2 - T^). (00 ^ 0.5: a = 0.75, 6 = 0.71 and 1 ^ 0.5: a = 0.70, 6 = 0.66. 



In order to see where the Kovacs-hump stems from in the present model, let us consider the 
populations in some more detail. In FigjTj we plot the populations pT2ie,ti;t) defined implicitly 
in eq.([8]). After a quench from Tq = 00 to Ti = 0.25 the temperature is changed to T2 = 0.5 
at the time ti needed for ET^{ti) = E^ given by E^ = — 2(T^. For t = 0, the population 
PT2{^,ti',t = 0) = ^1) is just the same one as immediately before the temperature change 

from Ti to T2. This distribution (the narrowest one in FigjT]) is slightly asymmetric and quite 

10° 
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Figure 7: Populations for the Kovacs protocol after the temperature jump from Ti to the final 
temperature T2 = 0.5 after t = 0, 0.01, 0.03, 0.06, 0.1, 10^ from bottom to top as indicated by the 
arrow. The dotted vertical line indicates the demarcation energy e£)(ti) = —Tiln^Kooti) and the 
symbols are for t = (black) and t = 10^ (red). 




narrow compared to the equilibrium distribution, which is in complete harmony with the results 
for a quench as discussed in ref.[28]. Note, that by definition of ti the distribution pTi(e,^i) yields 
E^ = —2(7^ for the mean energy. 

For finite times t a strong tailing of the distributions Pt2{^i ti] t) at higher energies is observed. 
Conversely, for low energies all curves but the one for the longest time superimpose indicating 
that relaxation towards the equilibrium distribution takes place only after very long times. This 
means that essentially all states with energies smaller than the so-called demarcation energy 
^niti) = — Ti In (^ooti) (shown as the dotted vertical line) are frozen for a long time|2H [30l [3T] . 
For intermediate times (all but the longest) the distributions qualitatively show the appearance 
of a superposition of two distinct distributions, one essentially given by pTi(e,^i) and another one 
with more weight at higher energies and increasing intensity for longer times. We mention that 
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this behavior is reminiscent of that found in the exponential trap model [TT]. 

Prom the behavior observed in FigjTj the emergence of a Kovacs-hump becomes clear. The 
mean value of the energy, -Et2(^i +t), will first increase and only after the low-energy states start 
to relax, it will decrease again and finally yield AE{t) = 0. 

The above discussion indicates that the population can be given as a superposition of two 
distributions. This means that the Kovacs hump can approximately be written as a differ- 
ence of two relaxation functions. In the following, we will show that for the trap model 
such an approximation can be deduced from the master equation, pT{e,t) = —KT{e)pT{e,t) + 
p{e)Jde'KT{^')PT{^','t), which is to be complemented by an initial condition ^^-(e, 0). In the 
second term, the gain-term, of the master equation, we approximate pT{€',t) by its long-time 
limit, the equilibrium distribution p^(e'). Furthermore, we utilize detailed balance in the form 
Jde' p{e)KT{e')p'^{e') = Jde' p{e')KT{^)PTi^) = i^t{^)Pt^{^) to obtain the approximative equation 
PT{e,t) = —KT{e)pT{^,t) + KT{e)p'^{e). The solution of this inhomogeneous equation poses no 
problem and for the Kovacs protocol, i.e. the initial condition pT{e,0) = pTi{e,ti), one finds 

p^P-(6,ti;t) = Nit,t,)-' [p^^ie) + {pT.(e,ti) - ^^^(e)} e-'^-.(^)*] (12) 

where N{t,ti) = 1 + Jde — p^(e)|e~''^2('=)* is a normalization constant. The dependence 

of N(t,ti) on ti can usually be neglected for a large part as can be seen from FigjTj For small 
energies one has PTii^^h) ^ Pi^i^) ^^^^ N(t,ti) ~ N{t) = 1 — /(iep^(e)e~'^^2('^)*. For 
e > e£,{ti), the relaxation is fast meaning that e~''^2('=)* ~ and consequently N{t,ti) ~ 1. 
Numerically, neglecting the dependence of N{t,ti) on ti introduces a maximum error (for times 
t ~ ti) on the order of 10%. For other temperatures Ti and T2, the neglect of the ti-dependence 
of N{t,ti) introduces an error of similar magnitude. 

For the distribution of the populations following a quench, the approximative solution has been 
compared with the exact numerical solution of the master equation in ref. [28] and good agreement 



has been found. We note, that a similar agreement is observed for the pT2{^,ti;t). Eq.(12) shows 
that the Kovacs hump approximately will be given as the difference of two functions. With the 
further assumption N{t, ti) ~ N{t), one finds that one function depends on ti and t and the other 
one only on t, 

AE(t)^0(ti;t)-^(t) (13) 
with 0(ti;t) = N{t)-^ Jde-e-pT,{e,ti)e-^T2^'^^ and 4j{t) = A^(t)- Vcie-e-p^^(e)e-''^2(^)*. A behavior 



of the hump similar to eq.(13) has been proposed earlier in ref.[TT] (eq.(27)) and it is also similar 



to what is obtained in linear response theory, cf. below. 



2. 'Inverse' Kovacs experiment 



Next, let us consider the Kovacs effect starting from a low temperature instead from a high 
temperature. The protocol then assumes the following form, see FigjSj As in the case of a quench 
from To > Ti, T2 (cf. Figs-3j|4), we first determine the time ti at which ET^{ti) = and the 



results are shown in FigjOj In contrast to the behavior for the standard Kovacs protocol, in this 
case ti changes monotonously as a function of temperature difference. This is because Ti is the 
highest temperature and therefore the relaxation time at Ti is the shortest. Only for Ti — )■ T2 it 
takes a long time to reach E^ because the energy has to relax almost to its equilibrium value. 
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Figure 8: Sketch of the protocol for the Kovacs effect starting from a low temperature. 
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Figure 9: ti versus (Ti — T2) (T2 = 0.5,0.6,0.7) for the Kovacs protocol starting from Tq 
(upper red lines) and Tq = 0.3 (lower black lines). 
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An example for the Kovacs-hump is shown in Fig 10 for Tq = 0.2 and T2 = 0.5. The behavior is 
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Figure 10: Kovacs-hump for Tq = 0.2cr and T2 = 0.5a and various values of Ti(> T2), Ti = 
0.55, 0.6, 0.7, 0.8 from top to bottom. 

- apart from the sign - similar to what is observed in Fig|5]for a quench from a high temperature. 
The minimum AEk increases with increasing IT2 — Ti| and tk decreases. 

For a direct comparison with the results of the standard protocol, in Fig, 11 we plot AEk and 
Tk versus IT2 — Ti\ for Tq = 0.8 and for Tq = 0.2. The most important feature of this plot is 
the similarity in the characteristics of the Kovacs hump for high and low initial temperature Tq. 
In both cases, one has power-law behavior of AEk and tk as a function of \T2 — Ti\ for small 
temperature differences. However, also the well known asymmetry of the relaxation [32] becomes 
obvious via the much longer time tk for the jump from a low temperature. 
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Figure 11: \AEk\ (upper panel) and tk (lower panel) vs. \T2 — Ti\ for Tq = 0.8 and Tq = 0.2. 



The fact that there is no qualitative difference between the two protocols seems somewhat 
puzzling on first sight. As mentioned already above, in case of an up-jump, the distribution 
of energies, pj-{e,t), shows a two-peak structure for intermediate times. The same behavior is 



observed when considering ^7^3(6, ti; t) for various times as shown in Fig, 12 From this plot one 
can see how the states with energies smaller than the demarcation energy euiti) relax as time 
progresses. The reason for the fact that this two-peak structure of the populations is not reflected 
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Figure 12: Populations pT2{^,ti;t) following the inverse Kovacs protocol for various times t = 
0,10^10^3 -10^10^3 -10^, 10^106 (from top to bottom). The intensity of the low-energy peak 
decreases with increasing time, indicating that states with e < en (dotted vertical line) relax. The 
symbols are for t = (black) and t = 10^ (red). 



in the shape and position of the Kovacs hump lies in the inherent averaging performed when the 
energy is calculated according to eq.(|6|. As in the 'standard' protocol, the Kovacs hump will 
approximately be given by a difference of two relaxation functions, cf. eq.(13). 

The expression given in eq.(13) along with the inspection of the populations for the two pro- 
tocols can be used to give a natural explanation for the sign of the Kovacs hump. For short 
times, the functions (p and ip can be written as 4>(ti]t) ~ — t /c/e ■ e ■ pTi(e, ii)'tT2(e) and 



m 



-C/T2 



t /(ie-e-p^(e)KT2(e)- Now, assuming a not too low temperature T2, one can further 
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approximate KT2(e) — t^oo^^'^^'^^ (1 + /92(e — E^) + ■■■). With this, one finds 



where the variances are defined as = Jde-{e — -p^i^e) and Vq = Jde-{ 
(remember that pT^(e,ti) = px^i^.ti^t = 0)). Thus, a glance at FigjT] and Fig, 12 shows that for 
the standard protocol one has AE{t) > and for the inverse protocol AE{t) < holds. 

Even though the Kovacs effect reflects the complexity of the relaxation properties of the system 
in some sense, its overall behavior does not depend very sensitive on the underlying distribution 
of relaxation times. This flnding might also explain why different models yield very similar results 
concerning the Kovacs effect. According to what has been discussed above the appearence of the 
Kovacs hump is a consequence of the fact that some energies have already relaxed at the time ti 
of the jump from Ti to T2 while others are still out of equilibrium. From FigjT] and Fig 12 and the 
discussion of the role of the demarcation energy it is clear that in both cases one might be able to 



model the Kovacs hump as the difference of two functions, cf. eq.(13). This results in a situation 



that is very reminiscent of what is observed in linear response theory [22] which we will discuss in 
the following. 



3. Linear response theory 

Using the superposition principle of linear response theory, one can calculate the response of the 
Kovacs protocol for |To — Ti| — )■ and IT2 — Ti| — )■ in the same way as in the calculation for a 
single temperature step, cf. eq.Q, with the result (cf. ref.[22]): 



Er.ih +t) = E't^ + [e'tI - E'^] Ceq(ti + t) 



^T2 



Ceq(t) 



(14) 



where Ceq(t) is evaluated at Tq and again AE{t) = -Et2(^i + 1) — E^. As mentioned already in 
the Introduction, one is concerned with the difference of two relaxation functions evaluated at the 
same temperature but at different times. Of course, linear response theory only holds for very 
small temperature jumps. Instead of showing results for the Kovacs hump here, we note that 
these qualitatively are very similar independent of the magnitude of the temperature differences. 
Only the quantitative features depend on the temperature regime considered. If one starts from a 



low temperature instead, the results are very similar. In Fig, 13 we show the results for the peak. 



lAE"^! (eq.(lO)), and the time of its appearance, tk (eq.(ll)), as a function of the temperature 
difference for some very small and some large temperature differences. From this plot it becomes 
apparent that the Kovacs effect does not change its qualitative behavior when leaving the linear 
response regime. The small differences in the values for tk for the positive and the negative 
difference show the limitations of the linear response regime and thus the onset of the asymmetry. 

Eq.([l4)) directly shows that a Kovacs hump will only be observed if the energy relaxation 
function Ceq(t) decays in a non-exponential fashion, as discussed in ref.[22]. Also in the nonlinear 
regime we flnd that the Kovacs hump can approximately be described as the difference between two 
relaxation functions, albeit they are not the same. Therefore, here a Kovacs hump can be observed 
also for exponentially decaying functions, if the relaxation times of (p and ip are different. However, 
for the trap model, all relevant functions decay non-exponentially. A notable exception from this 
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Figure 13: AEk, tk vs. IT2 — Ti| for different initial temperatures Tq 

behavior seems to be provided by tlie relaxation function in the p-spin model if one chooses 
p = 3, as is done in the quoted model calculations fT5]. In this case, the long-term relaxation in 
equilibrium is given by a single exponential function (Ceq(t) ~ e~*), but the relevant relaxation 
functions in the nonlinear regime are non-exponential. However, the relevance of such a behavior 
for the interpretation of the Kovacs effect in supercooled liquids is not obvious. Another model 
that shows exponential relaxation is the quoted two-state model[T9]. Here, the non-equilibrium 
situation is imposed by quenches with a finite cooling rate, a situation that we do not consider in 
the present paper. We expect that both models do not show a Kovacs effect in the linear response 
regime. 

If we relax the requirement that the temperature jump Tq — )• Ti be small, but consider only 
small IT2 — Til, the Kovacs protocol can be viewed as a linear response experiment performed 
on an aging system. In this case, one can use the formalism described in detail in ref. One 
just has to recognize that now the temperature jump from Ti to T2 plays the role of an applied 
external field, H = (T2 — Ti)/Ti (and the parameters determining the coupling of the field are to 
be chosen as /i = 1 and 7 = 0). The response is then determined by the energy autocorrelation 
function, where now this function is evaluated under the non-equilibrium situation. 

III. Conclusions 

In the present paper we have investigated the memory effect or Kovacs effect in a simple model 
for glassy relaxation, the Gaussian trap model. This model is interesting because the system 
can reach thermal equilibrium at all temperatures and it captures some aspects of the relaxation 
observed in real glass-forming liquids and computer models. 

The relaxation of the system following the Kovacs protocol shows a behavior that is similar to 
what has been found in earlier investigations on other models and also experimentally. However, 
we find one distinctive feature in our model calculations that has to our knowledge not been 
observed previously. The temperature dependence of the time ti, which is defined by the condition 
-E'ti(^i) = E^, shows a non-monotonic behavior and it first decreases, goes through a minimum 
and then increases again as a function of the temperature difference IT2 — Ti|. We suggest that 
IT2 — Ti| has to be rather large in order to observe this effect. Thus, we would propose to perform 
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experiments for larger temperature differences. However, a quantitative estimate of tlie required 
temperature differences is hampered by the fact that this depends on the relaxation strength and 
thus on the variable considered. 

For an up-jump and then an 'inverse' Kovacs protocol we find a behavior that is very similar 
to that following a quench, the main differences being a different sign of the Kovacs hump and a 
monotonic behavior of ti. The two-peak structure of the energy distribution does not affect the 
Kovacs hump significantly due the inherent averaging in the calculation of the mean energy. 

In both cases the different relaxation behavior for states with energies smaller and larger than 
the demarcation energy can be used to rationalize the origin of the Kovacs hump. Low energy 
states are still frozen at ti, the time of the temperature jump to T2, while large energy states 
respond quickly. From this fact it follows that the Kovacs hump can approximately be described 
as a difference of two relaxation functions, AE{t) ~ (j){ti,t) — i(j{t). The sign of the hump can 
be rationalized by an explicit short-time expansion. The concept of the demarcation energy is 
quite general and it can be used to understand many features observed in the aging dynamics 
of glass-forming liquids including the main features of the Kovacs effect. We interprete the fact 
that many models exhibiting a broad distribution of relaxation times qualitatively show a similar 
behavior as stemming from the distinction of slowly and fast relaxing entities. It is possible 
to relate the mobility of the entity to its potential energy. This establishes a direct connection 
between the microscopic system and the energy- related trap model [26j. The similarities found 
in different models can be understood as stemming from the fact that the averaging inherent in 
the evaluation of bulk properties is not very sensitive to the detailed shape of the underlying 
distributions. Consequently, the Kovacs protocol does not seem very well suited to distinguish 
among different models for dynamic heterogeneities in glass-forming liquids. 

When the temperature differences all are small one can use linear response theory for the 
calculation showing that the Kovacs hump has its origin in the subtraction of the energy autocor- 
relation function evaluated at two different times with different weights. We have shown that the 
qualitative behavior does not change for larger temperature differences. If only IT2 — Ti\ is small 
but |To — Til is arbitrary, one also can use hnear response theory but now one has to explicitly 
consider the non-equilibrium situation, resulting in a dependence of the energy autocorrelation 
function on two times. 

We conclude with noting that the Kovacs hump for a given quantity can be qualitatively 
understood in terms of the linear response behavior. In the nonlinear regime mainly quantitative 
differences are found but we do not observe any qualitatively new features. 
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